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ABSTRACT 

We report results of an extensive analysis of the X-ray nova XTE J1859+226 ob- 
served with BeppoSAX and RXTE during its 1999 outburst. We modelled the source 
spectrum with a multicolour blackbody-like feature plus the generic Comptonization 
model BMC which has the advantage of providing spectral description of the emitted- 
radiation properties without assumptions on the underlying physical process. The 
multicolour component is attributed to the geometrically thin accretion disk, while 
the Comptonization spectrum is claimed to originate in the innermost sub-Keplerian 
region of the system (transition layer). We find that XTE J1859+226 covers all the 
spectral states typical of black- hole sources during its evolution across the outburst, 
however during the very high state, when the disk contribution to the total luminosity 
is more than 70% and the root mean square variability < 5%, the high-energy photon 
index is closer to a hard state value (T ~ 1.8). The BMC normalization and photon 
index T well correlate with the radio emission, and we also observed a possible satu- 
ration effect of r at the brightest radio emission levels. A strong positive correlation 
was found between the fraction of Comptonized seed photons and the observed inte- 
grated root mean square variability, which strengthens the idea that most of the fast 
variability in these systems is provided by the innermost Compton cloud, which may 
be also identified as a jet. 

Key words: X-rays: binaries - radiative transfer- accretion, accretion discs - stars: 
individual: XTE J1859+226. 



1 INTRODUCTION 

X-ray novae, or soft X-ray transients (SXTs), are low- 
mass X-ray binaries which undergo a sudden and unpre- 
dictable increase in intensity up to 2-6 orders of mag- 
nitude in X-rays (|Liu et al.l 12001). These outbursts occur 
usually on time-scale of decades, even if there are cases 
like GX 339-4 where regular annual outbursts are seen 
l|McClintock fc RemillardfeuOtJ) . Their light curves are char- 
acterized by a typical rise-time of few days and an e-folding 
decay time of several tens of days. In these systems the 
accreting compact object is fou nd to be a neutron star 
(NS) or a black hole (BH) (e.g.. iTanaka fc Shibazaklll996l ; 
ICherepashchukll200Clk During the outbursts, thought to be 
originated by a dramatic change in the mass accretion rate, 
the X-ray spectrum evolves through different states, classi- 



fied according to the presence and relative prominence of 
a soft thermal component and a hard tail component (e.g. , 
I Zdziarski fc Gierliriskil|2004l; iRemillard fc McClintockll2006l: 
Icadolle Bel et al.ll20061 ; |Del Santo et al.ll2008l ; lLaurent et al l 



120111 ). The thermal component is commonly thought to 
be generated in the inner regions of the accretion disk, 
and the non thermal hard tail by Comptonization of 
the soft thermal photons b y a hot electron corona (e.g., 
iPoutanen fc Svensso 11 1996). The X-ray spectra of several 
SXTs may allegedly show a reflection component and/or 
emission lines and edges, due to reprocessing of the hard 
X-ray pho tons by the cold outer regio ns of the accretion 
disk (e.g. iMagdziarz fc Zdziarskil Il995l ). The power spec- 
tra derived from the X-ray light curve of SXTs in outburst 
are characterized by a red noise component with superim- 
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posed quasi-periodic oscillations (QPOs) at different and 
quickly varying frequencies. Black hole SXTs show generally 
no power above ~100 Hz and low frequency QPOs, while 
the power spectra of NS SXTs extend to higher frequen- 
cies and typically show kHz QPOs (e.g.. Ivan der Klis 2000; 
iBelloni et afll2002h . 

The X-ray outburst is accompanied by an optical 
outburst, with typical ratio between X-ray and optical 
luminosities of ~1000. The reprocessing of the X-ray 
emission from the outer regions of the disk is a pos- 
sible way to explain the UV, optical and IR emission 



vay to explain tne U V , optical and IK em ission 
van Paradiis fc McClintockl 1l995l; lEsin et al.l [2000; 



Coriat et a l. 2009; Gandhi et al. 201lT). By means of optical 



spectroscopy in quiescence, it has been possible to determine 
the mass function of several of these systems, which is the 
most direct way to establish if the compact accreting object 
is a NS or BH. Radio emission has been detected from sev- 
eral SXTs in outburst, and in some systems, often referre d 
to as " microquasars" , jets could be resolved |Fender| [2006). 
Simultaneous X-ray and radio observations of SXTs, albeit 
the global picture is not yet completely clear, seem to show 
some properties shared by almost all sources. In the canon- 
ical X-ray hard state, a steady and powerful jet is always 
observed with moderate velocity (P ^ 2), an associated flat 
radio spectrum (a ~ 0) extending beyond the radio band 
and linear polarization level of a few percent. During the 
source transition from hard to soft state, major radio flares 
are observed with ejection events stronger (P 2) than the 
hard state. The radio emission gets subsequently quenched 
as the source has reached the soft state, and newly turns-on 
in the soft to hard transition at later stages of the outburst 
l|Fender et al.ll2009l ). Albeit classical studies of jets in SFXTs 
have been performed in the radio band, in the recent years 
it has been shown that the contribution to the optical and 
infrared emission can be carried-out for a significant fraction 
by t he jet, both in blac k hole and neutron star X-ray bina- 
ries jRussell et alj2 006). Strong indirect evidence for optical 
and infrar ed fast-variab l e emi ssion from a jet has bee n re- 
ported bv lMalzac et al.1 (12004 ) and lHvnes et all i|200rj ). 

In this paper we report and discuss the results of an 
observational campaign of XTE J1859+226 performed with 
BeppoSAX and RXTE during the source 1999 outburst. In 
Section we give an overview of the source properties from 
the literature, in Section [3] we report information on the 
periods of observations with the two satellites and on the 
data reduction procedure, in Section [4.21 we show the main 
results of the spectral analysis, in Section [5] we discuss the 
implications for the source accretion geometry, some con- 
siderations on the relation between the high-energy spectral 
features and the radio emission, the main observational dif- 
ference with the systems hosting a neutron star and the 
origin of the fast temporal variability. Finally we draw our 
conclusions in Section [6] 



2 THE OUTBURST OF XTE J1859+226 

The transient source XTE J1859+226 was first detected by 
the All-Sky-Monitor (ASM) aboard of the Rossi X-Ray 
Timing Explorer (RXTE) satellite on 1999 October 9 (MJD 
51460) at the Galactic coordinates I = 54.05°, b = +8.61° 
with an intensity of 160 ± 15 mCrab in the 2-12 keV band 



l|Wood et alj|l999h . The ASM light curve (Fig. [TJ shows a 
structure reminiscent of a Fast-Rise and Exponential-Decay 
(FRED) shape typical of SXTs, with a rise-time (from 10% 
to 90% of the peak) of about 5 days, a maximum intensity 
of 1.4 Crab and a decay-time (e-fold) of ~ 23 days. This 
behaviour is however not smooth, and in the temporal inter- 
val from MJD 51468 to MJD 51479 during the decay phase, 
a strong flaring activity occurred. 

The outburst of XTE J 1859+226 was also revealed 
by the BATSE experiment onboard the Compton Gamma 
Ray Observatory up to 200 keV, with a peculiar be- 
haviour: the hard X-ray flux reached its maximum (MJD 
5146 2) while the soft X-ray flux (ASM) was still increas- 
ing (|McCollough fc WilsonTll999l '). Analysing the soft and 
hard X -ray temporal behaviour and comparing it with other 
SXTs, iBrocksopp et al.l (|2002| . hereafter B02), concluded 
that the hard state observed in this source at the outburst 
rise (despite the name of "soft" X-ray transients), rather 
then being a peculiar behaviour of XTE J 1859+226, can be 
an ubiquitous feature for this class of objects. 

The tempor al properties of the s ource were thoroughly 
investigated by ICasella et al.l l|2004l . hereafter C04) using 
extensive RXTE observations (from MJD 51462 to MJD 
51546). The Power Density Spectra (PSD) obtained with 
the Proportional Counter Array (PCA, 3-50 keV) showed 
three main types of low- frequency (1-10 Hz) QPOs, accord- 
ing to the centroid frequency value, the coherence parame- 
ter Q, the QPO amplitude, and the shape of the underlying 
noise (see Table 1 in C04). For C-type QPOs (the ones hav- 
ing Q 10), a strong anti-correlation was found between 
the frequency and the root mean square (RMS) variability, 
while the RMS amplitude of all the three classes increases 
with energy, flattening above 10 keV, ruling out a disk ori- 
gin. 

The source was also detected during its quiescent 
state in February 5, 2003, by the Chandra observa- 
tory; the 0.75-2.7 keV spectrum could be fitted by a 
simple power -law with a poorly constrained photon in- 
dex r ~2.4 (|Tomsick et al.l I2003T ). The optical counter- 
part of XTE J1859+226 was discovered on 1999 Octo- 
ber 12 at the position of R.A.(2000) = 18 h 58 m 41 s .58, 
Dec(200 0) = +22°39'29".4, w ith an 7?-band magnitude of 
~ 15.1 (|Garnavich et all 1 19991 ) and was observed with dif- 
ferent o ptical and infrared telescopes, including HST and 
UKIRT dHvnes et al.ll2002i ). The optical spectrum was char- 
acterized by a strong blue continuum wit h weak emission 
lines t ypical of an X - ray n o va in outburst llGarnavich et al.l 
1 19991; Wagner et al.1 1 19991 : ISanchez-Fernandez et al.l l200ll : 
Izurita et alj|2002h . Optical light curves present a FRED be- 
haviour similar to that observed in the soft X-ray band, with 
maximum magnitude V ~ 15 and quiescent magnitude V ~ 
23, a decay rate of 0.018 ± 0.002 mag per day. Evidences of 
3 mini-outbursts near q uiescence ( San chez- Fernandez et all 
l200ll : IZurita et al.|[2002n . a reflare at about 50 day after th e 
maximum and possible superhamps (jUemura et al.l |2004| ) , 
were also found. 

A strong radio counterpart of the source was discov- 
ered on 1999 October 11 at a flux densi ty of ~ 10 mjy at 15 
GHz with the Ryle Telescope and VLA l|Poolev fc Hiellmingl 
Il999h . The source position of R.A.(2000) = 18 h 58 m 42 s .0, 
Dec(2000) = +22°39'10".0 is consistent with the posi- 
tion of the optical counterpart. The radio light curve of 
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Table 1. Log of BeppoSAX TOO observations of XTE J1859+226. For each NFI, the left column reports the time exposure (ks), the 
right column the mean count rate (counts s~ 1 ) in the given energy band. 



TOO 


Epoch 


MJD 




LECS 




MECS 




HP 




PDS 








0.15-4 keV 


1.5-10 keV 


i 


i-30 keV 


15-200 keV 


1 


Oct. 15, 1999 


51466 


off 




14.1 


181.1 ± 0.1 


14.7 


63.6 ±0.2 


7.1 


42.6 ±0.1 


2 


Oct. 22, 1999 


51473 


off 




16.1 


159.1 ± 0.1 


17.2 


48.8 ±0.1 


7.7 


37.9 ±0.1 


3 


Oct. 28, 1999 


51479 


0.4 


163.4 ± 0.7 


4.4 


223.6 ±0.2 


4.0 


45.8 ±0.3 


1.8 


29.9 ±0.2 


1 


Nov. 7, 1999 


51490 


0.7 


139.5 ±0.5 


28.8 


135.9 ±0.1 


31.4 


9.6 ±0.1 


13.6 


8.2 ±0.1 


5 


Nov. 19, 1999 


51501 


4.5 


73.9 ±0.1 


30.1 


81.3 ±0.1 


32.3 


4.3 ±0.1 


12.8 


8.2 ±0.1 


6 


Mar. 25, 2000 


51628 


18.2 


0.80 ±0.01 


40.8 


0.97 ± 0.01 


45.2 


0.97 ± 0.10 


21.3 


1.14 ± 0.01 



XTE J1859+226 was characterized by five major radio flares 
superposed to a somewhat global fading behaviour (B02). To 
first approximation, the sequence of radio flares appeared to 
be correlated with the local minima in both the X-ray inte- 
ty i 

energies (jFender et alJl2009ft . 

The following parameters of the system were inferred 
from multiwavelength observations d uring outburst and 
quies c ence: orbital period of ^9.15 h l|Garnavich fc Quinnl 
|2000|; IFilippenko fc Chornockl 120011) high inclination i ^ 
60° l|Hvnes et all |2002|; IZurita et alj 120021), mass function 
f(M) = 7.4±l.lM m (IFilippenko fc Chornockll200ll). source 
dista nce of ~6-ll kpc I Hvnes et al.l |2002| ; IZurita et al.l 
12002ft . 

The values of the orbital period and mass function were 
recently update d to 6.58 ± 0.05 h and f(M) = 4.5 ± 0.6M Q , 
respectively, bv lCorral-Santana et al. M2OIII), based on pho- 
tometric and spectroscopic observations of the secondary 
star. The sour ce distance was instead updated to d = 
4.2 ±0.5 kpc bv lShaposhnikov fc Titarchukl (|2009l . hereafter 
ST09) using RXTE observations, and applying the scaling 
method to the F/QPO or T/disk normalization diagrams. 
With the same procedure, these authors also determined 
the BH mass to be M hh = 7.7 ± 1.3 M©. Fina lly, the mea- 
sured optical extinction E(B-V) = 0.58±0.07 (|Hvnes et all 
2002) is consistent with the ave rage N_g value of 2.21 xlO 21 
cm derived from radio maps (|Dickev fc Lockman|[l990l) . 



3 OBSERVATIONS AND DATA ANALYSIS 

Six Target of Opportunity (TOO) observations of 
XTE J1859±226 were carried out from October 15, 1999, 
to March 27, 2000, with the BeppoSAX /Narrow-Field In- 
struments (NFIs). They consist of a Low Energy Con- 
centrator Spectrometer (LECS, 0.1-10 keV, Parmar et al. 
1997 ), two Medium Energy Concentrator Spectrometer 
(MECS, 1.3-10 keV, Boella et al. 1997), a High Pres- 
sure Gas Scintillator Proportional Counter (HPGPC, 3-120 
keV, Manzo et al. 1997) and a Phoswich Detection Sys- 
tem (PDS, 15-300 keV, Frontera et al. 1997). During the 
first two TOOs the LECS was switched off. The reduc- 
tion of the LECS, MECS and HPGSPC data were per- 
formed with the standard FTOOLS/SAXDAS procedures 
whil e data from the PDS were redu ced with the XAS pack- 
age (|Chiappetti fc Dal Fiumdl 19971 ) . The energy bands used 
for spectral fitting were limited to those where the response 




Figure 1. RXTE/ ASM light curve of the entire XTE J1859±226 
outburst and zoom of the of the interval MJD 51460-51495 show- 
ing the flaring activity of the source before the beginning of the e- 
folding decay phase. The dashed vertical lines indicate the epochs 
of the BeppoSAX observations, while the two arrows show the in- 
terval of data analysed by ST09. 



functions are well known, i.e., 0.15-4 keV, 1.5-10 keV, 8-30 
keV and 15-200 keV, for the LECS, MECS, HPGSPC and 
PDS, respectively. In Table 1 we report the epoch of each 
TOO together with the exposure time of each instrument 
and the average count rates in the above mentioned energy 
bands. A systematic error of 1% was added in quadrature to 
the statistical uncertainties of the spectral data, on the ba- 
sis of the calibration results obtained with the Crab Nebula, 
that was observed on September 25, 1999, and on April 10, 
2000, i.e. 20 days before the first TOO and 15 days after the 
last TOO. We used the standard response matrices officially 
delivered for each of the four NFIs, except for the LECS data 
of T003, T004, and T005, as where the source count rate 
is higher than about 50 cts s _1 , specific response files must 
be created using the LEMAlQ routine. 

We accumulated total energy spectra for each TOO and 
for each NFI. This procedure was justified by the fact that 
the spectral analysis showed no systematic residuals between 
the data and the adopted models, meaning that the time 
scale of significant source spectral evolution was longer than 



http: / /www. asdc.asi.it/bepposax/software/saxdas/lemat. html 
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the instrument integration time. Normalization constants, 
allowed to vary in the recommended ranges, were introduced 
to take into account known differences in the absolute cross- 
calibration between the detectors 

A total of 131 RXTE observations from the public data 
archive were analyzed, from 1999 October 11 (MJD 51462) 
to 2000 July 24 (MJD 51749). The energy spectra and power 
spectral density (PSD) were extracted using the standard 
"RXTE Cookbook' recipes in HEAsoft 6.10 analysis software 
release. All RXTE data products were corrected for PCA 
deadtime. The background for energy spectra was produced 
by pcabackest tool version 3.8. The PCA response files were 
created using the latest release of pcarmf version 11.7 and 
energy-to-channel conversion table version e05v04 B The 
PSD was computed in the 0.01-64 Hz frequency range and 
normalized to give a RMS deviation per Hertz. In the fit- 
ting procedure, we ignored PCA channels corresponding to 
energies below 3 keV and above 50 keV, while the HEXTE 
spectra were constrained in 20-250 keV range. A free cross- 
calibration factor was used between PCA and HEXTE nor- 
malizations resulting in an expected value of ~ 0.85 for both 
HEXTE clusters. A systematic error of 0.5% was injected in 
the fitting process. 

We also checked how many single RXTE pointings oc- 
curred during each BeppoSAX TOO and found the follow- 
ing results: two for T002 (40124-01-22-00G and 40124- 
01-15-00), one for T004 (40124-01-44-00) and three for 
T005 (40124-01-51-00, 40122-01-04-00 and 40122-01-04-01, 
respectively). We subsequently tried to produce strictly si- 
multaneous BeppoSAX spectra in the time interval corre- 
sponding to the each single RXTE IDs. However, because of 
Earth occultation time, the effective exposure time for each 
NFI in these sub-intervals was too low, and eventually it was 
not possible to extract spectra with a statistics good enough 
to perform an accurate joint BeppoSAX / RXTE analysis. 

The spectral analysis was p erformed using the XSPEC 
package v. 12.5 (Arnaud 199g), and all the quoted errors 
for the spectral parameters correspond to 90% confidence 
level for a single parameter (A^ 2 = 2.71). When fitting the 
BeppoSAX spectra, alone or joined with RXTE, the inter- 
stellar absorption along the source direction was modelled 
using the elemental abundance by Anders & Ebihara (1982) 
and opacities from Morrison & McCammon (1983). On the 
contrary, because of the lack of data below 3 keV in the 
PCA, the spectral analysis of the RXTE spectra was per- 
formed by fixing the interstellar absorption to the galac- 
tic value Ah = 0-2 x 10 22 cm -2 along the source direction 
ijDickev fc Lockman|[l990l ). 



4 RESULTS 

4.1 Spectral modeling 

We modelled the source X-ray spectrum with a 
photoelectrically-absorbed two-component model, con- 
sisting of a soft and hard feature. For the soft component, 

2 ftp:/ /ftp. asdc.asi.it/pub/sax/doc/software-docs/saxabc- 
vl.2.ps.gz 

3 http: //www. universe. nasa.gov/xrays/programs/rxte/pca/ 
doc / rmf/ pcarmf- 11.7 
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Figure 2. Absorption-corrected deconvolved spectra in EF(E) 
units, best-fitting model DISKBB+BMC, and residuals between the 
data and the model in units of a for BeppoSAX TOOl and T004 
of XTE J1859+226. The single spectral components are also plot- 
ted: DBB {dotted line) and BMC (dotted-dasjed line). The emission 
line, marginally evident in the first and second BeppoSAX TOO, 
has not been included here. 



we tested both a simple blackbody ( BB), and the 
mult icolour-BB model diskbb in XSPEC l|Mitsuda et all 
1 1984 1 

The latter is obtained from the convolution of sev- 
eral BBs at different temperatures with Xbb oc i? _3//4 , 
where R is the radial d istance from the central object 
(Shak ura fc Sunvaevlll973r ). It is worth noting that the inner 
disk temperature (fcTm) and projected inner radius (Ri n ) of 
diskbb have to be considered as order of magnitude, rather 
than exact values, because of the intrinsic simplifications of 
the model which do not include colour-correction factors or 
physical prescriptions for the inner boundary condition (e.g., 
iGierlihski et alj 120081 ). Despite its simplicity, the diskbb 
model proves to be good in fitting the soft X-ray spectra of 
both NS and BH bin ary systems l|McClintock fc Remillardl 
120061 ; iLin et all [20071 ). 

The high-energy component of the source X-ray spec- 
trum was instead fitted with the BMC model whose emerging 
spectral shape is given by 

F(E) = ^ J [S(E)+AS(E)*G(E,E )], (1) 

where G(E,Eo) is the Green' s function (GF) of th e 
Comptonization energy operator (|Titarchuk et alj Il997l ). 
and S(E) is the input BB spectrum. The quantity A/ (A+ 1) 
in equation (fl} empirically describes the Comptonization 
fraction of the input seed photons. 

The free parameters of BMC are the seed BB tempera- 
ture fcTbb, the energy index a of the Comptonization GF, 
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Table 2. Best-fitting parameters of the model wabs*(dbb+bmc*highecut) for the six BeppoSAX observations of XTE J1859+226. The 
projected inner disk radius and luminosity are computed assuming a distance of 4.2 kpc. 



Parameter 


TOOl 


T002 


T003 


TOQ4 


TOQ5 


T006 



Nn (xlO 22 cm~ 2 
kT in (keV) 
ili n Vcos i (km) 
kT bh (keV) 

r 

log(A) 

^bmc 

E f (keV) 
^disk/ ^tot 
^0.1-200 keV 
-^0.1-200 keV 

X 2 /dof 



31+ 013 

061 +0.03 
u - D1 -0.03 

35+j 

1 00+ 04 
1 - uu -0.04 
n ic+0.04 
z - lo -0.04 

20+ ' 04 
U - ZU -0.03 

120+ 004 
u.izu_ 004 

77 +7 

"-6 

0.56 
3.67 
7.81 
131/109 



32+° 13 
u -° -0.13 

62+ 03 
32+| 

1 01+ 004 
J - ul -0.04 
n ic+0.04 
z - lo -0.04 

17+ - 04 

-0.03 

o-io±g;gg* 

0.56 
3.23 
6.90 
99/105 



+0.01 
0.01 
+0.03 
-0.03 
+2 
2 

+0.06 
0.05 
+0.08 
0.09 
+0.06 
0.05 
+0.01 
-0.01 
+57 
32 

0.64 
3.76 
8.02 
133/158 



0.28 
0.77 
24 
1.08 
2.17 
-0.13 
0.11 
130 



0.30 
0.71 
24 
0.89 
1 

-0.75 
0.047; 
159 



+0.01 
0.01 
+0.03 
-0.03 
+2 
2 

+0.04 
0.03 
+0.07 
-0.08 
+0.09 
0.07 
+0.008 
-0.008 
+72 
39 

0.77 



2.37 
5.05 
134/150 



0.27 
0.59 
28 
0.77 
1.81 
-1.04 
0.038 
> 130 
0.75 
1.59 
3.40 
160/146 



+0.01 
0.01 
+0.03 
-0.03 
+3 
3 

+0.04 
0.03 
+0.12 
0.13 
+0.09 
0.06 
+0.007 
0.008 



1.06 



+0.09 
0.08 
+ 0.03 
0.02 
+84 
33 
+ 0.04 
0.04 
+0.05 
0.05 
+0.14 
0.13 

+0.16 



0.39Z 

0.16j 
95^ 
0.32j 
1.87j 
0.19 1 * 



.0.11 



x 10" 



0.67 
0.10 
0.21 
162/138 



<Tn units of 10 -8 
''In units of 10 37 erg s 



cm 

1 



the Comptonization fraction log(A), and the normalization. 
The latter quantity, is defined as C n = L39/D1Q, where L39 
is the seed photon source luminosity in units of 10 39 ergs 
s _1 , and D10 is the distance in units of 10 kpc, respectively. 

The GF of the BMC model only depends on the spec- 
tral index a, namely G(E,Eo) oc E a+i for E < Eo, and 
G(E,E ) oc E~ a for E > E . As the GF is a broken-PL 
with no cut-off, the output spectrum at high energies is yet 
a PL. When the data show the presence of a high-energy 
rollover, this must be taken into account by multiplying BMC 
by an e-folding factor oc e~ E ^ Ef (highecut). 

In this sense, BMC is the most general available Comp- 
tonization model: it just measures the relative contribution 
of the direct and Comptonized component and the slope of 
the Comptonization spectrum, no matter which is the un- 
derlying physical process and thus having no limitation in 
its applicability. The origin of the observed features is de- 
manded in a second step to a scientific discussion. 



4.2 BeppoSAX and RXTE spectra 

The spectral analysis was first carried out with BeppoSAX , 
because of its broader energy band (0.1-200 keV) with re- 
spect to RXTE (3-200 keV), particularly in the critical re- 
gion where interstellar absorption is present. 

The BeppoSAX spectra of the source show the presence 
of a thermal component at low energies (Fig. [2]). However 
a simple blackbody (BB) model does not provide in general 
acceptable fits, while the diskbb better describes the data. 
Thus we adopted wabs*(dbb+bmc) as basic spectral model 
for the continuum emission. 

During TOOl, the presence of a high-energy cut-off 
in the spectrum was clearly evident. Indeed, the basic 
model does not provide acceptable fit (x 2 /dof= 690/110), 
while multiplying BMC by an e-folding factor, the model 
wabs*(dbb+bmc*highecut) gives a satisfactory result 
(x7 dof = 131/109). 

During T002, the spectral properties of the source were 
very similar to TOOl (Fig. [3]). In particular, the model 



a. 




Energy (keV) 

Figure 3. Spectral evolution of XTE J1859+22 across the whole 
outburst. The spectra of the six BeppoSAX TOOs and that cor- 
responding to the first RXTE observation (MJD 51462) are re- 
ported. 



wabs*(dbb+bmc) without cut-off gives x 2 /dof=451/160, 
while with cut-off we find x 2 / dof =99/105. 

The spectrum of the source during T003 at higher 
energies was slightly harder than that of TOOl and 
T002. In particular, the model wabs*(dbb+bmc) provides 
X 2 /dof=167/158, while wabs*(dbb+bmc*highecut) gives 
X 2 /dof=133/158. The PCI for the inclusion of the cut-off 
components is ~ 0.08, and it was obtained through the F- 
test routine for discriminating among two differ ent spectral 
mode ls for the case of multiplicative component (|Press et al.l 
Il992l ). 

During T004, the presence of the high-energy cut-off 
also proved to be marginal (Figs. [2] and |3}- The models 
wabs*(dbb+bmc) and wabs*(dbb+bmc*highecut) in- 
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deed give x 2 /dof=166/152 and x 2 /dof=135/151, respec- 
tively, with PCI=0.2 in both cases. 

The gradual shift towards higher energies of the spectral 
cut-off became evident at the time of T005, with the model 
wabs*(dbb+bmc) yielding x 2 /dof=164/147, and model 
wabs*(dbb+bmc*highecut) yielding x 2 /dof=160/146. 
Actually, it was only possible to put a lower limit on E{ > 
130 keV. 

At the time of the last BeppoSAX observation, T006, 
the source was moving towards quiescence with a signif- 
icant drop in the bolometric luminosity. At low energies, 
the statistics of the source was significantly worse with re- 
spect to the previous observations (Fig. [3]). At high ener- 
gies, the spectrum was dominated by a pure PL compo- 
nent and, unlike the case of T004 and T005, we could 
not put even a lower limit on the value of the cut-off energy 
E c . Thus, we just used the model wabs*(dbb+bmc) which 
yields x 2 /dof=162/133. 

During the first two TOOs, we found some evidence 
for the marginal presence of a Fe emission line (see Fig. [2]). 
Modelling the feature with a Gaussian component added 
to the continuum, the fit improves to x 2 /dof=lll/107 and 
X 2 /dof=85/103 for TOOl and T002, respectively. In both 
cases, the line width a was fixed to 0.2 keV, as leaving it 
free during the fitting procedure, the centroid of the line 
shifts towards too low values (E c ~ 6 keV). H aving in mind 
the issues reported bv lProtassov et al] (|2002T ) about the use 
of F-test in the case of emission lines, we empirically tested 
the significance of the feature. Specifically, we generated 10 4 
fake spectra for all the NFIs using for the continuum the 
best-fitting model wabs(diskbb+bmc*highecut) without 
the Gaussian emission line, and with the same time exposure 
of the true spectra for each NFL 

We then fitted the simulated spectra first with the 
continuum model alone, and later with inclusion of the 
line (with a = 0.2 keV). For each simulated spec- 
trum, we computed the value F12 = (x? — Xa)/(yi — 
V2)v2 /xh where J12 is the prob ability distribution func- 
tion |Bevington fe Robinson|[2003l ') . Then, we computed the 
number of spectra N satisfying F12 > F°2 B , where F°£ B was 
computed on the true data. We found that the fraction of 
simulated spectra obeying the above condition is 38 x 10~ 4 
for TOOl and 140 x 10 -4 for T002, respectively. 

These values represent the empirically-estimated prob- 
ability of chance improvement (PCI) for the addition of the 
emission line. 

As far as the RXTE spectra are concerned, we note that 
ST09 considered only the PCA data in their spectral anal- 
ysis of the early part of the outburst of XTE J1859+226. 
The lack of spectral coverage below 3 keV and the lower 
energy resolution of the PCA with respect to the LECS 
and MECS experiments onboard BeppoSAX , allowed ST09 
to fit the source spectrum with a simple photoelectrically- 
absorbed GAUSSIAN+BMC model. However, the BeppoSAX 
low-energy data clearly showed the presence of a soft ther- 
mal component, which actually carries-out a significant part 
of the source luminosity (see Table 0). Thus, driven by the 
BeppoSAX results, we tested both wabs(bb + bmc) and 
wabs(diskbb + bmc) models also on the RXTE spectra. 
We found that in both cases, they satisfactorily fit the data, 
strongly reducing the significance of the Gaussian compo- 
nent. For consistency with the BeppoSAX results, we thus 



have chosen as reference model wabs(diskbb + bmc), fixing 
the photoelectric absorption to the value across the source 
direction Ah = 0.2 x 10 22 cm" 2 . However, given that more 
than one hundreds RXTE observations were performed dur- 
ing the source outburst, instead of reporting all the best- 
fitting parameters in a table, we decided to plot them as a 
function of the MJD overplotted with the BeppoSAX results 
(Fig.©. 



5 DISCUSSION 

5.1 Spectral evolution across the outburst 

The extended temporal coverage of XTE J1859+226 by 
RXTE together with the BeppoSAX observations allowed 
us to deeply investigate the spectral evolution of the source 
during the whole outburst (Fig. [31). 

The X-ray continuum is characterized by the presence 
of a soft thermal component, described by a multicolour BB 
model (diskbb) and, at higher energies, by the convolution 
of a seed BB spectrum with a Compton up-scattering GF 
(bmc). Even if the data did not allow to fix equal BB and 
diskbb temperatures (fcTbb and kT[ n , respectively) in the 
fit, we interpret them as having a common origin in the 
accretion disk. In particular, as fcTbb > kTi n , it is worth at- 
tributing the BB temperature to the innermost part of the 
accretion disk. We refer to this region, where the gas gets 
hotter and most of Comptonizat ion takes places, as the tran - 
sition layer or Compton cloud |Titarchuk fe Fioritd 12004'). 
A similar spectra l mod elling using BMC was also adopted by 
iMontanari et all (|2009l . hereafter M09) in their BeppoSAX 
analysis of XTE J1650-500, even if in that case a simple BB 
could adequately fit the soft excess. 

In Fig. [3] we report the behaviour of the spectral param- 
eters of the best-fitting model wabs*(dbb+bmc*highecut) 
as a function of time (MJD), while in Fig. [5] we show the 
integrated (0.1-64 Hz) RMS and the soft (2-20 keV) and 
hard (20-200 keV) deconvolved X-ray fluxes obtained in the 
period where also radio observations were available (see Sec- 
tion E3J). 

As already reported by B02, XTE J1859+226 entered 
its outburst in a hard state, with the peak of the soft X-ray 
light curve delayed by a few days (Fig. [S[ see also Fig. 1 
in B02). At the peak of the hard X-ray luminosity (MJD 
51462) the source spectrum was characterized by a domi- 
nating thermal Comptonization (TC) feature, with cut-off 
energy below 100 keV (Fig. O and photon index F ~ 1.5. 
In this phase the disk contribution to the total luminosity 
was less than 20%, and the PSD showed high temporal vari- 
ability (RMS ~ 30%, see Fig. Q), both features typi cal of a 
canonical hard state |Remillard fe McClintockll2006l ). 

During the rising phase of the soft X-ray flux (Fig. 
[TJ), the Comptonization spectral index T progressively in- 
creased, reaching its maximum (~ 2.5) in coincidence with 
the soft X-ray peak. In the same period, the Comptonization 
fraction log(A) decreased, as well as the BMC seed photon 
temperature. On the other hand an increase of the inner disk 
temperature kT[ n occurred, together with a decrease of the 
projected inner disk radius. The simultaneous and mono- 
tonic decrease of the RMS is also indicative of a spectral 
transition from hard to soft state. 
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This trend of the spectral parameters is usually 
expected in the scenario of the truncated disk model 
ijDone et al.ll2007h . 

The subsequently observed flaring phase of the ASM 
light curve is claimed to mainly originate in the accretion 
disk, which dominates the spectral emission below 10 keV 
(Fig.EJ. Indeed between MJD 51466 and MJD 51480, fluc- 
tuations in the inner disk temperature (however not in its 
normalization) occur, and they are approximatively coinci- 
dent with the secondary peaks of the ASM light curve. 

Not surprisingly, any variation of the disk properties 
is expected to have physical consequences in the radiative 
processes of the Comptonization region: hardening and soft- 
ening of the high-energy X-ray spectrum are visible in Fig. 
[4] through fluctuations of V. Once again, changes of the ge- 
ometrical extension of the corona, in response to this effect, 
are still visible looking at the behaviour of log(A). Yet, dur- 
ing the flaring phase, significantly less visible changes are 
observed in the feTbb value or BMC normalization. 

It is of particular interest the source behaviour in the 
period between BeppoSAX T004 (MJD 51490) and T005 
(MJD 51501), during the smooth ASM exponential decay 
phase. In this phase, the disk contributed more than 70% to 
the total luminosity (see Figs. [2] and [3]), and, at the same 
time, RMS was < 5%, typical of a soft state. 

However, the photon index was V ~ 1.8, which is actu- 
ally a value closer to a hard state (see, e.g., XTE J 1650-500 
in M09). 

Note also that this time period also corresponds to the 
minimum value achieved by the Comptonization fraction 
parameter log(A), as shown in Fig. [4] thus presumably in 
this phase the Compton cloud reaches its minimum size. As 
also shown by the BATSE light curve reported in B02 and 
by BeppoSAX and RXTE in Fig. [5] while moving towards 
quiescence XTE J 1859+226 experienced an increase of the 
hard X-ray flux (> 20 keV) around MJD 51520, with a cor- 
responding increase of the RMS (~ 10%). 

This feature however does not necessarily implies a tem- 
porary soft to hard transition at all. Actually, Fig. U shows 
that the rise of the hard X-ray luminosity was accompanied 
by an increase of the Comptonization fraction log(A) and 
a steepening of the photon index T. At high energies, the 
spectrum thus became slightly softer but globally brighter. 
This trend appears puzzling when looking at the behaviour 
of the two soft seed photon components. Indeed, a slight in- 
crease of the disk temperature kT[ n occured a few days later 
(around MJD 51530) after the new rise of the hard X-ray 
luminosity, and appears not tightly correlated with changes 
in the BMC parameters log(A) and F. 

The only appreciable correlated variation is found in 
the BMC normalization (but not in feTbb), which actually 
also increased around MJD 51520 more or less tracing the 
behaviour of V. Thus, we interpret the rebrightening of the 
hard X-ray luminosity as an enhancement of the amount of 
photon supply from the innermost region (bmc normaliza- 
tion) such that this effect compensated, in terms of Comp- 
tonization energy flux, the steepening of the the photon in- 
dex r. 

Note that this is a different behaviour with respect to 
what happened at the beginning of the outburst, where the 
coronal behaviour appeared mostly dictated by changes in 
the diskbb component (see above). 



The increased time intervals between each RXTE point- 
ing for MJD > 51540 and the reduced statistical quality of 
the high-energy RXTE spectra, did not allow a thorough 
mapping of the evolution of the source spectral parameters. 

However, more than five months after the outburst on- 
set (MJD 51629), BeppoSAX observed the source when the 
0.1-200 keV flux was ~ 10" 9 erg cm" 2 s" 1 . The broad-band 
coverage and sensitivity of BeppoSAX allowed the unique 
opportunity to catch the source close to quiescence with an 
excellent data quality (Fig. [3]). In the phase down to qui- 
escence, the spectrum of XTE J 1859+226 is characterized 
by a strong PL component with T ~ 1.9 and no evidence 
of cut-off. The accretion disk contributes yet to a signifi- 
cant fraction of the total luminosity (about 67%), but now 
is much cooler (fcT bb ~ 0.3 keV, kT in ~ 0.2 keV), and the 
Comptonization fraction is at the same level of the rising 
phase of the outburst. 

The very low temperature of the direct disk component 
does not allow its detection by the PCA, and this explains 
the gap of data points in the two lower panels of Fig. [4] 
The spectral shape of the last BeppoSAX observation can 
be interpreted as a new significant increase of the Compton 
cloud characteristic size, together with a very high electron 
temperature, consequences of a strongly reduced efficiency 
in Compton cooling by the seed photons. 

5.2 Justification of the adopted spectral model for 
describing the high-energy emission 

The high-energy emission of BH sources in different spectral 
states has been fitted wit h several XSPEC Comp tonization 



models such as COM ptt dCadolle Bel et al.ll2006l). com pps 
l|Frontera et al.ll200lh , or eqpair (iDel Santo et aTTl2008T ) . All 



the above mentioned models contain a detailed treatment of 
the radiative transfer physics, but each of them present its 
own drawbacks. 

In particular , it is worth e mphasizing that both COMPTT 
l|Titarchuklll994h and COMPPS (|Poutanen fc Svenssor]|l996T) 
are static thermal Comptonization models with different 
ranges of applicability for the plasma electron temperature 
and optical depth (preferably low kT B and high r for COMPTT 
and the opp osite for compps). On the other hand eqpair 
(|Coppilll999T ) allows to compute Comptonization spectra for 
a plasma which can be purely thermal or hybrid (thermal 
plus non-thermal electron distribution). 

Thus, when one of such models is used to fit the source 
spectra at high energies, it is implicitly assumed that the 
spectral changes occur in the framework of static-only ther- 
mal Comptonization processes (comptt, compps) or due 
to a switch from static thermal to hybrid Comptonization 
(eqpair) . 

By fitting the RXTE spectrum of XTE J1859+226 dur- 
ing the bright hard state at the beginning of the outburst 
(MJD 51462, see Fig. with diskbb+compps model, as- 
suming spherical geometry of the Compton cloud, we found 
kT c ~ 30 keV and r ~ 1 (with Comptoniztion parame- 
ter Y ~ 0.2). Following the source evolution from the hard 
to the soft state we note however that, e.g., the fit of the 
BeppoSAX spectra provided an electron temperature kT e 
increasing from ~ 50 keV (TOOl) to ~ 200 keV (T005) 
with a corresponding decrease of the optical depth from r ~ 
1 to ~ 0.2, respectively, and Y m 0.3 — 0.4. 
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Figure 4. Evolution of the parameters of BMC and DISKBB models over the outburst of XTE J1859+226, for the RXTE {green points) 
and BeppoSAX (blue points) observations. The vertical dashed lines indicate the times of the five radio flares identified by B02 (see also 
Fig.O 



The x values were acceptable in all cases, but in our 
opinion it is difficult to explain such a behaviour during 
a hard to soft transition. In particular, the enhanced disk 
contribution to the total luminosity (see Table[5} is expected 
to provide strong Compton cooling of the coronal electrons 
and their temperature should decrease rather than increase. 

It is also worth mentioning that because of the RXTE 
and BeppoSAX high-energy threshold around 200 keV it 
was not possible to find evidence of any additional PL- 
like component above the energy rollover suc h as reported, 
e.g. in GX 334-4 by iDel Santo etU] l|2008l ) with INTE- 
GRAL/SPI. Thus, hybrid Comptonization (eqpair), rather 
than excluded, is not required by the data. The progressive 
increase of the high-energy cut-off during the hard to soft 
transition, together with the spectral index steepening, can 
be more easily explained by the transition from a TC to bulk 
Comptonization (BC) dominated state (e.g., M09). 



Unfortunately, specific models for BC in BH sources 
are no t yet available. The COMPTB model of iFarinelli et all 
(2008) is based on an assumption of sub-relativistic bulk mo- 
tion, expected to hold in NS systems due to strong radiation 
pressure from the NS surface, but not in BHs. 

Actually, our choice of BMC was mostly motivated by 
the fact that with its parametrization, it represents the best 
compromise between a phenomenological and physical de- 
scription of the high-energy source spectral behaviour, possi- 
bly avoiding contradictory results which may obtained, e.g., 
using static TC models. 



5.3 Correlation between X-ray and radio emission 

A thorough study of the correlation between the X-ray and 
radio emission of XTE J 1859+226 was previously carried- 
out by B02. In Fig. [5] we report the cumulative light curve 
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Figure 5. From top to bottom: radio light curve at all wave- 
lengths (see Section lOt . integrated 0.1-64 Hz RMS, 2-20 keV and 
20-200 keV fluxes in units of 10 — 9 erg cm -2 s _1 obtained from the 
deconvolved best fit model DISKBB+BMC using both BeppoSAX 
and RXTE. The time period of the plotted data corresponds to 
that of the available radio observations. 
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Figure 6. Spectral photon index V (upper panel) and BMC nor- 
malization (lower panel) as a function of the radio flux normalized 
at 15 GHz for XTE J1859+226. Data have been averaged over one 
day. 



of the source obtained from several observing facilities at 
frequencies ranging from 1.22 GHz to 22.5 GHz in the period 
from 1999 October 13 (MJD 51464) to 1999 December 13 
(MJD 51525), and we address the reader to B02 for details. 

The radio light curve shows a rather complex behaviour, 
with fluxes ranging from 0.5 to 100 mjy, and five major flares 
identified by B02 (drawn as vertical lines in Fig. [4] and [5} 
occurring over an otherwise globally fading behaviour. 

Starting from the results of the X-ray spectral analysis, 
we searched for correlated patterns between the radio emis- 
sion and the source spectral parameters in spite of integrated 
X-ray fluxes. 

To avoid biases in the obtained results, due to differ- 
ent observing times at different frequencies, we preliminary 
normalized the radio fluxes to those at 15 GHz, using the 
relation F(15 GHz) = F(v)(15 GHz/V) a mjy, with the spec- 
tral index a reported in Table 1 of B02. We then considered 
the normalized radio light curve and single parameters of 
the diskbb+bmc model averaged over half and one day, re- 
spectively, and we calculated the Spearman correlation co- 
efficient r. Because the number of points was more than 20, 
we also calculated the significance of the correlation using 
the associated t-variable of the Student distribution. 

The strongest X-radio correlation (> 99% confidence 
level) is found, over one-day time scale, between the radio 
flux and both the spectral photon index F and BMC normal- 
ization (iVbmc), with r ~ 0.6 in both cases (Fig. [6]) . 

On the other hand, the hard X-ray flux (given by the 
convolution term in equation [I]) is proportional to AWc, 
thus in turn its positive correlation with the radio emission is 
eventually expected as well. We find indeed that it is similar 
to that observed with the aim of BATSE data (see Fig. 8 in 
B02). 

It is interesting to observe that brighter and steeper 
X-ray spectra are associated with brighter radio events, with 
a possible saturation effect of the photon index F for increas- 
ing radio fluxes. 

As reported in Fig. [7] the same saturation effect is in- 
stead not observed in a diagram of V vs. AWc, a typical 
feature of most, albeit not all, BHCs (ST09). 

The index saturation in a F vs. AWc or F vs. ^qpo 
diagram is interpreted as an indirect evidence for the pres- 
ence of a converging flow (bulk motion) close to the cen- 
tral object in BH sources, an d it is expected to occur 
when the optical depth r ^> 1 (jLaurent fc Titarchukll2007l : 
iTitarchuk fc Seifinalr2009l ) . 

The thorough monitoring with RXTE during the whole 
outburst of XTE J 1859+226 (Fig. |4} seems to exclude the 
possibility of an incomplete sampling problem. 

As t oc M, one possibility, albeit speculative, is that in 
this source the accretion rate M does not exceeds the critical 
value above which the Comptonization parameter Y, and in 
turn the spectral index F, becomes insensitive to r. 

The apparent index saturation in the F vs. radio dia- 
gram is thus challenging in this context. 

Correlation patterns between the radio emission and 
the integrated R MS variability were investigated by 
iFender et all (2009) who found that the five radio flares iden- 
tified by B02 occurred around periods of reduced RMS. This 
is almost evident looking at Fig. [5] albeit the correlation is 
much less significant (r < 0.2). Also, the different behaviour 
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Figure 7. Spectral photon index V versus BMC normalization for 
the RXTE data set analyzed by ST09 (orange), plus the second 
part of the outburst (black) and the six BeppoSAX TOOs (blue). 
For clarity purposes we plotted only the points for which the 
ratio between the best-fitting value and the lc error determined 
by XSPEC is > 3. 
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Figure 8. Fraction of Comptonized seed photons, parametrized 
through log(A) in BMC, versus integrated 0.1-64 Hz RMS variabil- 
ity using RXTE data. 



of GX 339-4 and possibly XTE J1550-564 seems to point 
against a direct casual link. 



5.4 Transition layer as a formation site for the 
fast X-ray variability 

One of the most established observational results of BH sys- 
tems is that their PSD strongly correlates with the source 
spectral state. During the hard state, the RMS variability 
is high (up to 40%) and the PSD is described by a bro- 
ken PL with a flat plateau up to the break frequency. When 
the source moves progressively towards the intermediate and 
soft state, the RMS variability decreases below 10% and the 



PSD is a PL with index 1.0-1.5 (|Remillard fe McClintockl 
2006). The rapid variability down to ~ 10 2 s is probably 
mostly dictated by variations of the Comptonization region 
(TL), as the characteristic viscous ti me-scale of the accretio n 
disk is several hundreds of seconds (IChurazov et aljfeooil ). 

Models for explaining the rapid variability as a func- 
tio n of the source spectr al st ate have been propose d, e.g., 
bv lTitarchuk et al.l (120071 ) and llngram fc Pond (120111 ) in the 
framework of propagating perturbations in bounded config- 
urations. 

In order to strenghten the hypothesis that rapid vari- 
ability originates in the TL, in Fig. [8] we present a plot of 
the log(A) parameter of BMC versus the integrated 0.1-64 
Hz RMS of XTE J 1859+226. A tight positive correlation 
is observed between the two quantities for RMS > 5%. The 
Spearman rank correlation coefficient in this case is r — 0.95, 
while including the whole data it is r = 0.83. As it can be 
seen from the definition of the BMC model (Eq. [I]), log(A) 
quantifies the contribution of the Comptonized spectrum 
with respect to the input seed photon population. To first 
approximation, log(A) is an indirect measurement of the size 
of the Compton cloud. The fact that when log(A) increases, 
the corresponding source RMS increases as well, is thus a 
strong additional probe that it is the Compton cloud the 
site of origin of the rapid variabilty. 

On the other hand, unlike the case of the broad-band 
noise component, the origin of low-frequency QPOs both 
in NS (10-60 Hz) and BH (0.1-10 Hz) sources, is still 
controversials, as either m odels proposing a Lense-Thirring 



disk precession origin (e . g., Stella fc Vie tri 1998; S tella et al.l 



ik precession origin (e .g 

ll999l ; ISchnittmanll2005l ; ISchnittman et al.1 [20061 ) or magne- 
toacoustic oscillations o f the Comptonization region (e.g., 
iTitarchuk fc Osherovichl 1 19991 ; I Wagoner etHI l200ll ) have 
proved to be in good agreement with observations. 

A tight correlation between the spectral photon index 
r of the Compton cloud and characteristic C-type QPOs, 
which have variable frequency (0.1-10 Hz) and high coher- 
ence parameter (Q > 10), has been observed in a large 
sample of BH sources by ST09. Additionally, C04 (see their 
Fig. 3) found a tight anti-correlation between the RMS and 
^qpo for QPOs of the same type. These observational fea- 
tures, together with the fact that usually QPOs show greater 
amplitude at higher energies, apparently would favour the 
Compton cloud, rather than the accretion disk, as the region 
where QPOs form. 

However, it is worth noting that the above mentioned 
correlations are not in contradiction with a tilted disk origin, 
if most of the disk thermal soft radiation is intercepted and 
Comptonized by the hotter inner region. This interaction 
indeed, would shift to higher energy the soft photons mod- 
ulated at the Lense-Thirring disk precessional frequency, at 
the same time preserving the F vs. ^qpo and RMS vs. ^qpo 
relations. 

It is also worth mentioning that no correlation with the 
RMS was observed by C04 for A-type and B-type QPOs (7 
< v < 8 Hz, Q <3 and 4.5 < v < 6.5 Hz, Q > 6, respectively) 
and that the F vs. QPO correlations diagrams for BH mass 
determination with the scaling method were produced by 
ST09 using only Type-C QPOs. 

This pattern behaviour could be indicative of a different 
physical origin for the various types of QPOs with Type-C 
ones associated to magneto-acoustic oscillations of the TL 
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l|Titarchuk fe Fioritol [20041 ) and Type-A and Type-B ones 
originating in formation site either away from the TL or in 
some physical process (e.g., jet/corona connection) almost 
independent on the geometrical variations of the TL itself. 



6 CONCLUSIONS 

The joint BeppoSAX / RXTE analysis that we performed 
on the X-ray Nova XTE J1859+226 during its 1999 out- 
burst has demonstrated the importance of getting simulta- 
neous spectral and temporal informations during the source 
evolution. The extensive observational campaign performed 
by the two satellites indeed allowed us to map thoroughly 
XTE J 1859+226 from the onset of its outburst down a phase 
close to quiescence. We modelled the hard X-ray component 
with the BMC model which has the advantage of being a fully 
generic Comptonization model, suitable to be applied in any 
physical regime because the shape of the up-scattering GF 
function in BMC can be considered as general as possible. 

The source experienced all the canonical spectral states 
of BHCs from hard to soft (Fig. [3J). During the soft state, 
with a disk contribution to the total luminosity higher than 
70% and an RMS < 10%, the high-energy spectral index, 
with r ~ 1.8, was actually closer to typical values of the 
hard state. 

A significant correlation was observed between F and 
the radio emission, with a possible effect of index satura- 
tion at the brightest radio emission levels. Additionally, we 
found a tight correlation between the Comptonization frac- 
tion of the seed photons, parametrized through log(A), and 
the source RMS, which allowed us to strengthen the claim 
that most of the source variability originates in the tran- 
sition layer, where also Comptonization is the dominating 
process. This correlation is new and can be of key interest 
in understanding the accretion geometry of X-ray binary 
systems systems hosting either a BH or a NS. 

We thus strongly suggest to test it on a large sample of 
sources, which is possible due to the large amount of archival 
RXTE data. 
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